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ABSTRACT 

We construct compactly supported wavelet bases satisfying homogeneous 
boundary conditions on the interval (0.1). The maximum features of multires- 
olution analysis on the line are retained, including polynomial approximation 
and tree algorithms. The case of //^ ( (0, 1)) is detailed, and numerical values, 
required for the impleraentation.are provided for the Neumann and Dirichlei 
boundary conditions. 
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INTRODUCTION 

Wavelet bases are often presented as a powerful too! to perform the ap- 
proximation and the numerical iosohitikm of partial differential equations 
Indeed, thanks to rero moment and localisation properties, wavelet spaces 
are self adapted tc the solution and therefore may allow fast and accurate 
resolution. In the last few years different algorithms have been successfully 
tested on linear and non linear equations [9- [12]. Nevertheless, most of the 
addressed problems where posed in the periodic framework which circum- 
vents the difficulties generated by general boundary conditions, but which 
makes nearly impossible the treatments of real problem is. 

It s known (3) that under very general hypotheses, one can only consider 
homogeneous conditions. Therefore, we are driven to the construction of ap 
proximatkm spaces for homogeneous functional spaces on the interval. For 
example, for problem involving homogeneous Dirichlet conditions on [0,1], 
i.e ti(0) = «(1) ~ 0, the solution u could be reached in the Sobolev space 
^([0. !]) = (u € f J ([0, lj). u(0) ~ tt( 1 ) =0} and this leads to the construc- 
tion of approximation spaces for /f^([0. 1]). Various const roc lions have been 
prop*, ed (see Auscher (2), in [4] and references in it) but are not, according to 
their authors, numerically tractable. The proposed eonstructi n follows the 
one of Cohen et al. ((lj). After a short recall of this preliminary construction 
we introduce and ftaaJvxe tire homogeneous space construction. Numerical 
details are provided 


I THE PRELIMINARY CONSTRUCTION 
OF COHEN, DAUBECHIES AND VIAL 


This construction of wavelets on the interval, ([!’}. is derived from the com 
partly supported wavelet multi resolution analysis on the line introduced by 
t D&ubechies |5]. 

In the case of JL J I JR) this multi resolution analysis is classically given bv 
a sequence of closed subspares V, satisfying: 

i) •• C V-i C Vo r V\ C v 3 c . = {0} and U,es V ■ — 

L l {B) 


\ 


Moreover, e&cl V) is spumed by the translations of the dilated version of a 
fixed function, the scaling function i.e V, * — k), k € X}. 

Here, the family {^^(y. — k). k € X] is orthooormal and <p is on the 
one hand compactly supported and on the otlwr hand such that its Fourier 
ransform <?> satisfies the Strang and Fit approximation rules of order N l 

PI- 

^{2kr) = 0. k€X\{ 0}, n = 0 N~ 1. (1) 

One consequence of (1) is that the family, {&. ~ k) t k f_ iT), can 
reproduce locally the polynomial of degree at most N - 1. 

The support of d is the interval [—A' -f 1, A r j and the regularity of o is 
asymptotically C° xV [5]. Moreover. 6 is solution of the following scaling 
eg nation: 


<jHx) - y2 h^>( 2.r — k). (2) 

Jt* ~ ft -f l 

The detail spaces H’ ; are defined as the orthogonal complements of V, in 
Vi+>i *-e, 

if, = v, t , n (v;) i (:t) 

and. thanks to i) 

»€* 

The essential feature of multi resolution analysis (see Y. Meyer [ill) is that 
3c such that Vj € X 

W, = spaniV 11 ^’.- kl keX). 

Again, the family {2 ,/2 a(2 J . — k). k £ 2t] is orthonormal. The function 
V' is here a compactly supported wavelet and is obtained from the folkswing 
dttail equation 


,Y 

T(r)-- ^ gn^i'Ir k) 

*»■ A * I 


(41 



Moreover suppi^) = supp{<t>) and has the same regularity as o. In addition, 
because of the approximation properties of Vq and the definition of H 0 , p 
has got .V vanishing moments, ie: 

fx'p(: r)=0 

Finally, the family {<£;.*(*) = 2 lf7 p(2*x — &)} is an unconditional basis for 
various functional spaces such as Holder spaces C*(H) or Sobolev spaces 

H'{iRy. 

In [1], the goal of I. Daubechies et al. was to construct a family of wavelet 
basis on the interval [0.1] able to characterise I 3 ([Q, 1]). //*f[0, 1]) or C*([0. 1)) 
while preserving the most attractive properties of muhi '-solution analysis 
of L*( JR), despite the lack of shift invariance of i*((0, 1]) (this is not the case 
for the con struct ion* of P. Auscher ( ^ M)). 

We give in the following paragraph an outline of the construction but the 
reader should refer to [lj for details. 

The construction is performed in two steps as follows. 

The first step consists in defining suitable subspaces of £ 2 ([0, 1]) from 
a basis essentially constructed from the translated version.' >f a rescaled 
function, while the second step consists in the construction of the detail 
spaces with the same requirement. 

More precisely, in the first step, V'([0, lj) is instructed as follows: 

Thanks to the compact support of for large enough values of j and 
k — .\\ .,2* - .V — 1, the support of the functions - k) is included in 
[0,1]. Therefore, tb corresponding functions may be used a* the interior 

functions of Vy((0, lj) and the set ♦/ -= {+(‘2 3 x- k), k — A', ^ - A 7 — 1 } (/ 

stands for interior) is then defined. To fully define V'((0, lj). A ccf gt functions 
are added at each boundary of [0,1] to complete the basis ♦/. These two 
families of A functions, ♦k.o = { 9®*, k = 0, ..., A-lj and *e , = k = 

2-’ — A 2 } — 1} are constructed to have minimal support and, such that 

the order of approximation, (A), related to the interior functions is kept. 
In other words, all polynomials of degree lent than N - 1 should be locally 
expandable as a linear combination of the basis functions of V,([0, 1]). Let us 

’We remind that * € iR, f bekmgs to H‘{T) if aad only if C,,f * i/WiV 4 n J )* < 
that foi 0 < <» < 1 , / € C° if and only if | f(r > c h) -f{x)\ < C|h|* tor every 
x .h >a K the contt-ant C not depending o»i jr and h 


,1 


recAll the construction of ♦$,<) (the same applies to *e,\). The edge functions 
o are defined for k ~ 0,..., .V — 1 as the restriction to [0,1] of a specific 
linear combination of the family {^(a- - ft)such that 0 € supper *)}. 
More precisely. 


i.V-2 

“ y /7 ]£ C) * + n ~ X + i) X|o,i] (5) 

n»* 

The supports supp{^ 0 } k ) are staggered, i.e, = [0. ) and <p®* 

is a polynomial of degree k on the interval [O.^j (see Figure 1 for example). 

Indeed. n k $(x - n) is a polynomial of order k [1] and Q) is a polynomial 
in n of degree k. 

By construction, $£o X ♦/ but $g,o is not an orthonorraal family. An 
orthouormaliration procedure using the Grain- Schmidt algorithm is then pet 
formed. Starting from v-i down to one obtains A' orthonormal edge 
functions k = 0, .... A' — 1} w-ith staggered support [0, .V -f 1] and still, 

is a polynomial of degree k. 

Finally. i j([0. 1]) is then by definition generated by the orthooormal family 
♦jg.o U $r U $k,; : 


V’.{ it). 1]) — span < 


{#*. k ~ Q , ..., V - 1} 
U 


k = - A' - 1} 

u 

( {*»*. y - 1} J 

with Ojk * 2-’ /2 P(2-'r — F), One gets 

v»((0,l]) c V»+i([0.1j)C - C V'(io.lj) C ...£*((0,lj) 


(6) 


where ; 0 is choeen so that supporting#) n ,upporfi$ t , ) = 9. As 0 ; *. the 
edge functions satisfy a modified scoltny equation (2). one writes 

.V-l A' + 2* 




C H + £ A fc.n^j41.a * * 0, --.A - 1 

tiwC tic .V 


(7) 


The numerical value* of the coefficients n - 0, ...,A' 4 2i; Jb — 

0 A*— 1 } and (M.,,- n r-- 3*2* - S’~2k~X .. .2 J ~ 1; t * 2 >-,Y. . „2'~l} 
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for the right edge, are computed in (lj. 


The second step of the construction is the definition of a suitable basis 
for the usual wavelet space Wj([0.1j) = lj+t((QJ]) n (l;([0, l])) 1 . Thanks 
to (3) and to the compact support of %\ for large values of j , the family 
= {v-(2-’r — fc), k ~ - iV - 1} belong to Wj((0, 1]). Wavelets 

of this family constitute a first part of the basis of Wj([G, 1]) and are called 
the interior wavelets. Since Am(Hj([0. 1)) * 2-’, ,V other wavelets at each 
edge should be added to ♦j Again, we only recall the construction at the 
edge x — 0. Tis complementary wavelets are deduced from the definition of 
W0([O,l))as 


v-l 


V 0 ?* 


v j+r.r 






n*0 


r j ,n 


fc = 0, A' — 1 


( 8 ) 


where (.,.) stands for the scalar product of I*({G. 1]). By construction they 
are orthogonal to Vj([0, lj) and to the interior wavdets. Their supports are no 
longer staggered, but an iterative process described in [1] reduces the support 
of Vji to instead of [O,*—- 1 ]. The last step of this roastrucuon 

consists again of a Gram-Schmidt orthonormalmtion. Starting from k ~ 0 
up to A — 1 one gets A" orthonormalued wavelets for the left edge, { , k -- 

0, .... A — 1}. These wavelets are known through the coefficients {$ n = 
0. .... -V 4 2 k\ k — 0, .... A' - 1} that occur in the modified details eqvatton: 

A- 1 iV+2* 

V,.k = X. + X) k — 0,..., X — 1. (91 

H',((0. 1]) is therefore entirely characterised by 


Since 


w,((0, lj) - span 



k = 0...., A r - 1} 

1 


u 


W',i, k ■ 

= A', - A’ - 

1} 


U 


(vj>> k 

- 2' - X . , 2 J - 

1} . 


i 3 ((o.ij)« i r ((o,i!)(;D w-((o.i;i 

,1 Zfv 


( 10 ) 
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one gets an orthonormal basis of I J ([0, I]) as: 


t*0 M ~ 1} 


t-0 A’-l) \ 

u 


u 


U 

t«.V 2 J - A ; - 1} 

U 

J>*> 

U 

k {*»*< ** 2*-A\...,2»-l} . 


fe = 2 J - A\...,2 J - 1} J 


( 11 ) 


Remarks: 

As we have said before, these wavelets bases are very attractive because 
they preserve the main features of the whole line construction. More pres 
cisely, since the edge functions are finite linear combinations of some shifts 
of they have the same regularity. From their definition the edge scaling 
functions generate all the polynomials up to degree N ~ 1 which ensure an 
order <V approximation over all the interval, and the existence of S van- 
ishing moments for the edge wavelets. With these oscillations and enough 
regularity, these wavelets bam form an unconditional basis for the Holder 
spaces C*([0, 1]) [1]. The fast wavelet transform [10] which is essential for 
most numerical applications is preserved even near the boundary thanks to 
the modified scaling, (7) and detail, (9) relations. 


Our aim is to adapt this construction to obtain wavelet families generating 
functi nal spaces with homogeneous boundary conditions. More precisely we 
will consider the following constraints 

f <CLo) { 0) = / <rt,) (l) « 0, 

where /*'* is the i-th derivative of /. 

As will be shown, most of the above construction will be preserved as 
well as numerical efficiency. 


f> 



II MULTI RESOLUTION ANALYSIS WITH 
HOMOGENEOUS BOUNDARY CONDI- 
TIONS 


This section is devoted to the construction and the properties of compactly 
supported wavelet satisfying homogeneous conditions of tvpe 0) vs 

/< C£) >( 1)*0. 

The starting point has been described in the previous section and, keeping 
the same notations, we now assume that the compactly supported wavelets oo 
the line satisfy 0 < CLQ.CLl < N - 1 and r > max(CT0. CXl}. Therefore, 
the spaces l;{[0. lj) defined in (6) are included in C*((0, 1]) with r > a > 
m«x(CjL0, C£l). 

II. 1 Construction 

As in the previous section, we only focus on the left edge i~0 
According to (6), every function f } € V' ([0, 1]) is written 

A-l JJ-.V-l V - l 

fA*) * H W?,* + H H Wj.k- U?) 

k*: V k*V -«Y 

Moreover, only the lef edge functions ■#£ o are non zero around x — 0 and 
then 

/T>) = eV>(^.) < c “’( 0). (13) 

i*0 

Therefore one way to impose /' CLo) (0) = 0 is to enforce that all tl>e left edge 
scaling functions satisfy this condition. Following P. Auscher (in (2] ). this 
constraint can be related to a polynomial behavior 

Indeed, from the last section we learned that <fP k is a polynomial of de- 
gree A' - 1 on the interval (0. £], say for example p,.*(x) — a° k + a' k x -f 
- + a* k The CIO-nth derivative of at 0 is then equal to off®. 

Therefore (s^j*) (0) = 0 <=* a^j (® 0. The construction of scaling 

functions satisfying *®^(0) »» 0 is then equivalent to the construction of 
edge functions such that their restriction.* to [0.‘t~ 3 j as no component on the 


monomial x CLo . 


The first step of our algorithm is then to construct a family of jV edge 
scaling functions $*.o = k * 0, A — 1} and one of N edge vavelets 

^£,0 ~ k * 0 , N — 1} with the particularity that only one scaling 

function and one wavelet contain r Ci0 on their polynomial part. The second 
step is to remove the scaling function containing x CLo and to modify the 
corresponding wavelets. For simplicity we work on the interval [0,-f oo[ with 
a zero dilation scale (j = 0), omitted in the next notations. Moreover, we 
call p*(x) = a* + ajar 4- ... + the restriction of ^l(x) on [0,1]. 

We start with the first N edge scaling functions (5) of section I. They are 
defined with the coefficients c*° „ so that 

k SH-2k-2 

*?k(*) = £*L^(2x) + £ k ~ 0,.... N - 1 (14) 

»«C , Km S’ 


(see [1] for the computation of these coefficients). The following proposition 
tells us how to modify y* to eliminate x CL0 in the polynomials pt-k ^ CI O. 
We call *he new functions and Pk{x) ~ a\x' 

Proposition II. I The. family k = 0. .... A’ — 1} defined by: 



a cio,* “ HmCto+i 




is tvch «f i0 = Q,Vfr ^ CL0. 


S 



Proof: 

The existence of A* is always ensured for k f CIO since a° k k « 2~ k 
[1]. Because p* is a polynomial of degree F, there is nothing to change for 
k < CIO, and therefore 9° « ^ as well as p*(x) = p*(j t) for i < CIO. Given 
fc > CIO, let us suppose that Vf < kj ^ CIO, of 1,0 = 0. From relation (14) 
we obtain the following scaling relation for 9°: 

*=»(*) * -*><*«.,) #?(2*) 

+Ak > cL« V’cin(^ r ) a *,ct<Hi 9 cio«-i( 2*H ••• + a k.* 9»(2x) 

( 16 ) 

where 

k - 1 

Ak,cto — Q t.ctc ~ (°CLO,CLO ~ Q t,k) + 

Since <$(2x - n)|j 0 .t/j] s 0 far n > A\ the contribution of the third RHS 
term of 16 to is 0. Moreover, for 0 < i < k, 9?(2x) = pf(2x) are polynomial 
with no component on x CW) . Therefore, the contribution of p£ to i cio is 
entirely due to Ai,,clo and At,cLo — 0 is the condition we are looking for, 
that completes the proof. ■ 

For k > CIO the supports of <f>\ are no longer staggered but. in compen- 
sation Vi, 0 < k < N - 1.9k!(c,»j is still a polynomial of degree k. Therefore 
the functions 9®.0 < k < A r — 1 are independent. Moreover, they are or- 
thonormal to the |4(x - n), n > A'} since they are linear combinations of 
the {9®, k * 0, ..., A’ — 1). 

Following the previous section we now orthogoualiae the family {9^.0 < 
k < A’ — 1 ), keeping the “monomial independence.* The only thing to do is: 
to exchange the piece of ■Pclo an ^ 9*-i before starting f the algorithm from 
index 0 up to .V-J. The result is an orthonormal family of N edge scaling 
functions {9° J ", k = 0, ..., .V— 1} with the particularity that only pf t _ l . the 

restriction of 9^ti °° [0.1), contains x CLo . They satisfy a modified scaling 
equation: 

V-J sv-j 

9?- x (x) - V H?^ m x (‘2z) + £ H 0 k n 4(2T - n) k = 0 - !( 171 

H»0 R».V 
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with = 0 for k - 0,..., A r - 2. 

The construction of the N edge seeding functions for the right edge comes 
from the seme Algorithm for the half line joo,0j. The (f* 1 . it = 2 J -.V + 
1. .... 2 ; •- 1 } ere inde [indent of x CLl and only comains this monomial 
on [0,1]. 

After a dilatation of 2 3 for the 0 and 1 edge functions and adding the 
2 J - 2A : interior scaling functions dy,*. one gets therefore a new orthonormal 
basis of Vy ( [0, lj), the space defined by (6). In this family, only $£10 (*«P- 
contributes to x cu (resp. x Cil ) on [0, 2~ J ] (resp. (1 - 2~ J , 1]). 

To perform our first step of construction we now continue by isolating a 
single wavelet containing x cu> on [0, 1/2}. 

As in the previous section N wavelets at each boundary should be added 
to the interior family ♦/ Focusing again on the left edge, we construct a 
first family following (8) as 


,v-i 


^k(*) = 9k' X (2*) ~ 53 (^k 1 ( 2jr )^n A (*)) '•fn V) * ~ 0 A’ - MB) 


\*Q 


Again, each function vJ>*(r) is polynomial on the interval (0,1/2]. 

However, since for all it, depends on 9^-t,, all the v s ° contains the 
monomial x CLC and are therefore not suitable for our first step (we remind 
that we want to construct a family of edge wavelets such that only one 
contains x CLo on (0,1/2.). Still, from (17) and (18) we deduce a modified 
detail equation for these functions that writes 

= 53 l (2x) + 53 #,d(2x - n) 0 ,A’ - 1. (19) 

#1*0 g»iV 


1 he following proposition tell us how to transform the functions i'® to reach 
our first step. 

Proposition U.2 T he family if E .o - (V*®. A- = 0 A’ - 1} gxvtn by: 


/ Vk - C ° ~ (HV v_i k = 0, .. 

I V % * V'fv-i 

is such that only thr rrstrietton of c 


A f - 2 

with m ■= 1 


^N-l.A-1 
,V_j to [0.1/2] contains . 


( 20 ) 
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Proof: In (19) the monomial i CL0 is present only in Writing the 

details equation for and canceling the coefficient of fives the desired 
result, a 


As previously, we apply a Gram Schmidt orthonormalisation that pre- 
se^-ea the above property. Indeed, starting from v® ** get an orthoootmai 
wavelets family #£ >0 = { V*® 4 ', k = 0, ..., jV - 1 } for which only , contains 
x CLi on [0,1/2]. These wavelet* are defined using the details equation 


!»“(*) * L - ») t - 0 ;V - 1(21) 

wK) »s.V 

It only remains to make this construction again for the right edge with 
monomial x CLl and to expand all the boundary wavelets of a factor 2^. To- 
gether with interior wavelets familv they form an orthonormal basis of 

W,(f0.1]). 

We have nov rei d**d our first step since we have constructed a basis 
of scaling functions for and a basis of wavelets for IV. (j0. 1]) such 

that in each family, only one function has a component on z CLo on [0.1/2] 
and only one function has a component on i CLl on [1-1/2]. 


As announced, we now perform the second step of our const ruction by 
removing the function on the left edge and the corresponding ones, 

'jpJ.'i.y, for the right edge. 

The last technical point is the modification one should make to the wavelet 
space. We have the following proposition: 

Proposition II. 8 Define the sukspact V,([0. lj) as 

v:,([0. 1]) a l'([Q. 1]) - *pan{$X^i v^i_ lV } (22) 

Replace the fw rot rlcts v®‘y-i «nd m the families o and V Bl by 



„ -o.j. , . ; o,j 

a v,..v-i + "Vjji-i 





( 23 ) 
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with o and b solutions of 


| <i J + 6* =« 1 


(24) 


«n<f £ ’ Wk<jor3 of the some set of equations with coefficients for the right 
edge. 


Then the new family »•* «* orthonormal basis «/$")( |0, 1)), 

the orthogonal complement o/^((C, 1]) in V>n({0. lj). Moreover every scaling 
function o/Vj(jO, l)) and every wavelet o/U' ; ([0, lj) satisfies the homogeneous 
boundary conditions /* CM *( 0) = ) « 0. 


Proof: We prove only the result for the left edge. 

Let us first recall that s?° v..j and aw respectively two bah* func- 
tions of V,([0.1]) and U'([6,l)) v and that V)([0J]) JL H ; ((0,lj). ©? is 

then orthogonal t-o all the other basis functions of Vj([Q, lj) and Wj([Q, 1]|. 
Moreover,! |0ylj|i = 1 if and only if a*-f fr 3 =■ 1. The same argument holds for 
©1 and therefore, with the new definition of and ♦£,», oU ♦ / U » 
is a family of 2 J ortbonocmal functions. 

Using the scaling (17) and detail (21) equations we get 




JA’-J 




daO 


mi rnS 


Taking into account (24 1. we get that that 0® is independent of rViA-i and 
consequently belongs to V^ +1 ([0, lj) Since the orthononnal collection ♦g.o'J 
$ / U generates a closed subspace of Vj+i((0, 1)). orthogonal to f ; ((0, lj) 
and of dimension V = d»mVj + |((0, 1]) — lj), it is by definition 

Hj([0. lj) the or t honor mal complement of 5^((0, 1]) in Vj 4l (!0, lj). 

« 


Remarks: 

All these edge functions have the same regularity as the initia] scaling 
function d. 
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Moreover the 2 A - 2 edges wavelet* constructed before we remove ^'y., 
and have conserved their A' vanishing moment*. The modified 

wavelets 0® and 0] belong to VV’ ; ((0. 1)) and are therefore orthonormal to 
all the polynomials included in V,((0, 1]). Bu' they have no reason to verify 
(0®,* rw ) — (0],x Cil ) s* 0. since these monomials haw been excluded from 
the edges of V>([0, 1]). Hence only one vanishing moment for one wavelet at 
each boundary h r -s been lcet. 

As in the mitial construction, the modified scaling and detail relations 
insure that fast algorithms related to the different basis projections are avail- 
able. 

At this point however, we don’t know exactly what kind of space the 
multi resolution family Vj((0. lj) *oproxirna*e. This is the purpose of the 
next subsection. 


II.2 Approximation results 

We now check the intuitive result that ’he wavelet basis derived from the 
last construction is an orthonormal basis for suitable homogeneous spaces 
on [0, l].We give a complete proof for the Dirichkt homogeneous boundary 
conditions /(0) ~ /(l) = 0. i.e, CZ.0 - CL] — 0, corresponding to Hq ( [0, lj). 

Using our construction for CLO — CLl — 0. we first obtain a subspace 
\) ({0, 1]) defined by the octhom i nal basis #f >0 U ♦; U $£. i. with the partic- 
ularity that only two scaling functions are non zero at the edges. It is known 
that under some specific conditions (see the previous section! V'([0. 1 ) be 
longs to a multi resolution analysis of the Sobolev space // k ( } 0, 1]). 

Let vis simplify the notations and write — ,p® 4 and Then 

we have the following result: 

Proposition II. 4 Let V,([0. 1]) Ik the sub apace spanned by the orthonormal 
basis 


M.0i — i) LK^ n- -v } LKv?;,n-v } 

Assume these scaling functions hate enough regularity, to tnro/rr 


U^do.'Oi-H'do.ii, 

)>JJ 


n 


and that only and art non zero at 0 and 1 Then 

U 'i((0- 1!) = UmdoaD-^ni^..,;^,.,,.}) = HU 10. lj) 

J>JC )>» 


Proof: Take ft function / in //£([(), lj), H j ( / ) *ad n ; (/) its or- 

thononnal projection onto Vj((0. lj) said V,((0, 1]). We have Vo establish the 
relation 

Inn ||/-li ; t/)i!«. =0 (25) 

where the //‘-norm is taken as Following the 

density of Vj((0, 1]) in //‘((O, 1]) this is equivalent to 

lim !in,(/) - n^/HU: = 0 


Now using the ortho normal basis of Vj(fO, lj) and V^((0. lj) we haw 

since the support of the 0 and 1 edges scaling functions d- not overlap. 
Because — 1 and belongs to H l {[ 0. 1)'» (due to the regular- 

ity of the initial function £), we have 


litt/) - n,(/)ll»- < c,y (it/.y-J*..,); + K/.vU-.v'l) 

where Cj is a constant independent of j. Therefore we haw to check that 

lim y«/.el»-,)l - .Hm a’K/.vU-xH-o 

;-*+5C * j— 4-oc J 


To see this, we use the inequality 


^!as^A-i)!<^ll/i!H‘ (26) 

with Cj independent of j. which will be justified at the end of the proof. Let 
us take now a sequence of functions (/*)» € jv convergent to / with respect to 
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the //‘-norm, for example /„(jr) = f{x) ,_ij. Applying the list inequality 
to f — /, we get 

v\(f Call/ -All*,. 

It only remains to note that there exist an integer J, dependent of n, for 
which 

Vj > j, *\{f - /.,*?*.«>• “ W-v-’*-,)! 

Indeed, for fixed n we take J ruui that the supports of y>® iV _, and /„ do not 
overlap. 

Making n tends to -foe, and consequently j, leads to 

lira 2M</,<&v-,} i = 0. 

Obviously the same arguments holds for the scalar product 2-’!{/, v>“ 3J _ v )| 
and the proposition is proved. ■ 

We still have to establish the inequality (26). An integration by parts 
implies that 

where is a primitive of p-® s Since belongs to I J ([0, 1)) we deduce 

iiH?!|y < C\T-\ 

This assumption and the definition of the H x - norm leads to the desired result. 
Thanks to this proposition and the definition of W ^([0. 1’) (see prop II. 3) we 
deduce a decomposition of Hq([Q, 1]) in term of wavelet basis, 

«it(0. ill - 1»([0. 1|) 0 wyu 1’,). (27) 

)>x 

Remarks: 

The proof lor Neumann homogeneous conditions is similar and involves 
the // 3 -norr'i. More regularity is therefore needed for the basis functions and 
a double integration by parts to the inequality corresponding to (26). In that 
case, the approximated space is the strict subspace of f/ 5 ([0. 1)) defined as 

{/ e //' i (io,i)),/° ) (o) = / :i \n - o}. 


Some mixed homogeneous boundary cooditioi s, for example /( 0) = /( 1) — 
p l) \0) s /01(l) = 0. could also be addressed with a similar construction. In 
this case, two scaling functions at each edge are removed from V' ( [0, 1]) and 
are employed to modify the wavelets of H^([0. 1|). This construction leads 
to a characterisation of the functional space 

£/2((0,l])* {/€ /f 3 (!0.1)),/(0) = /(l) =0 and /^(O) « /™(1) - 0}. 
Since the left and right basis functions do not interact at scale j, different 
conditions could also be taken at 0 and 1. 

The following section is related to the numerical estimates related to our 
construction and to various topics connected to its application for partial 
differential equation problems. 

Ill NUMERICAL ESTIMATES 


This section is devoted to the numerical estimates related to our construction 
for two cases of homogeneous boundary conditions, i.e. the Dirichlet condi- 
tions and the Neumann > onditions. All the following computations have 
been carried out oeginntng with the initial compactly supported function d 
closest to linear phase constructed by I. Daubechies (6] with A ~ 4. Since 
no explicit analytic expressions exist, this function is defined through the 
filter coefficients h* used in the scaling equation (2). These coefficients are 
provided in [6]: 

= -.07576571478950. /»_, = -.2963552764600, h. t ~ .4976186676328 
ho = .8037387518051, h , = .29785779560531, k , = -.0992195435766 
/» 3 = -.01260396726203. = .03222310060405. 


The corresponding interior wavelet V’ is defined using the coefficients g n of 
the details equation (4} with = (— 


III. 1 Dirichlet Boundary conditions 

The application of the last section algorithm with CZ0 — CLl - 0 (Dirich- 
iet condition) leads to a multi I'esolution analysis of H,}([0. 1]) Three scaling 
functions and four wavelets have to lie added at each l»oundary (see section 
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II). These seeing functions Are solutions of * modified scaling equation (17) 
and are therefore characterized by the coefficients //£ H and //]! „. The cor- 
responding numerical estimates (computed on a 16 decimal digits computer 
with an error smaller than 10 ' u ) are listed in Table 1. The coefficients G° i n 
and G] „ which occur in the modified detail equation (2i) are listed in Ta- 
ble 2 and define completely the edges wavelets. All the following figures are 
obtained using the cascade algorithm [6]. The three x — 0 edge scaling func- 
tions. as well as the three r — l edge scaling functions are represented on 
Figure 2 at scale > * 0. The corresponding wavelets are plotted on Figure 3. 
Noth e that due to the lade of symmetry of the initial scaling functions and 
wavelets, the x ~ 1 edge functions can not be deduced from the x - 0 edge 
functions using a simple transformation. 


111. 2 Neumann Boundary conditions 

The same numerical estimates corresponding to th« Neumann conditions, 
i.e CL0=CL1-1, are listed in Tables 3 and 4. The Figures 4 and 5 repre- 
sent respectively the scaling functions and wavelets of this multi resolution 
analysis. 

Remarks: 

Some zero coefficients are provided in Tables 2, 3 and 4. They axe ex- 
pected as foil- . >: for instance in Table 3 N$ A *= H$ 2 = 0: since the scaling 
functions y^o.o for the Neumann conditions, is by definition constant on tb* 
interval [0.1], it does not dej>end on , and r ^ 2 which are respectively poly- 
nomials of order 1 and 2 on [0,1/2]: this leads H% { - //J 2 = 0 in Table 3 
Others zeros are expected using the same arguments. 

111. 3 Quadrature formula 

In order to use these wavelets basis for numerical purposes one question 
needs still to be answered. Given a function /, how can we define a projection 
Vj([0, 1]), i.e. how can we estimate a set of coefficients e,_* occuring in relation 
(12) and corresponding to /? The solution proposed here aims to compute 
an approximation of the orthogonal projection of / on lj([0, 1]) defining 
quadrature formula to estimate the coefficients fjj k m /.MU We define 
below a quadrature formula of order JY - 1 in the same philosophy as G 
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Beylkin et al. ([8]) or W. Sweldens et ad. ({13]). We we therefore looking for 
weight coefficients uy* such that 

*' ./(«,) (28) 

' <m0 

where the {a,, i = 0 A' - 1) are A' given points taken in the support of 

* and such that the approximation is exact for the polynomials of degree 
less than or equal to A' - 1. 

It appears that the weight coefficients uy* are the solution of the following 
linear system: 


/ * E “U*)' 1 = ' N ~ 1 (») 

J isO 

Hence we need to evaluate the first A r moments of every edge scaling function. 

Multiplying the modified scaling equation (17) by x 1 leads to the A r -- 1 
equations: 


sA’-i 

L 


,V - 3 S.V-J . 

/ xVSjx) - £ H?, / xVo.(2x) 4 £ H?„ / x'o(2i - n> * = 0 A' - 2 

J nmO J J 


(30) 


Since the moments of order / of the interior funct ion oi, / ar^<to’r ). can be esti- 
mated using the classical recurrence relation given in [8], (30) finally leads to 
the following linear svsiem .4 A? = h, where the X dimensional vectors A) and 

X,(i) - /xV'S^r) 

6; are defined as and 

m = rtlv ! *f../xV<2x-n).l- = 0 A-2 

and where the entries of the matrix A depend only on the „. We easily 
checked that this matrix is always nonsingular, (to see this , use the fact that 
E£? !«?.«! < 1 »>>«• 1 Vo * II t* = 1- 

We first provide the numerical values of the moments of order f.O < / < 3 
for A' = 4: 

M 0 = 1 OOOOOOOOOOOetOO M x = 1 1 934 5240«- (T2 

A/, = 2.1 1 177120S 1 8e-04 M* * 4.34510522842e-02 
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Then, the entries of A 7 for / = 0 . ...,A' - 1 for the 0 and 1 edge scaling 
functions corresponding to Figure 1 (Dirichlet boundary conditions) are listed 
in Table 5 and 6 . 

Using the values of these moments and the A' given points a,, we find the 
weights <a> for every edge scaling function solving the linear Vandermonde 
system (29). 

Remark: 

A quadrature formula of same order has to be used to estimate the in- 
terior scaling coefficients c,,* = / S < k < - .V - 1 to preserve a 

constant order of accuracy all over the interval. 


Acknowledgments: 

This work has largely profited from discussions with Philippe Tchamitchian 
and special thanks are due to him. 

IV CONCLUSION 

Compactly supported wavelets satisfying homogeneous boundary conditions 
on [0,1] have been constructed. All the tools required for the use of these 
functions for numerical approximation of partial differential problems have 
been detailed. 

Even if all this construction extends by tensor product arguments to 
higher dimensions, efficient handling of general open sets with boundary 
conditions is still an open problem. 
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LIST OF SYMBOLS USED 


The I^TpX code of the mathematical symbols used is given to clarify their 
identity: 

^([0.1]) : H‘l_0([0,i]) and a similar code for L a ff?) ; // 1 ([0, 1]1 ; 

i a ([0,l]);C'([(U]); 

V}([0, lj): \tilda{V>.j ([0,1]) and same for ^([0,1)) : \) : IV’j ; Vj([0. 1]) 

(J V', : \ovwrlin«{\bigcnp_{ j \in \2) V.j} and same for f) V}, {J H' ; 

d ( ’' ) (2*sr) : {\hat \phi}‘{(jt)}(2k\pi) ; 

4,*: \phi_{j,k}; 

$E.o : \phi_<E,0> ; ♦; : 

• \ t i lde{ \ vaxphi } ' {0 , \p«rp} _{ j ,k> ; The same expressions axe us- 

ing substituting \warphi by \pai and \Phi by \Psi. 

Q) : \8igl(*{n}_{k>\Bigl) ; (,): Mangle , \raaglw ; $ 

\«»ptys«t 

h k n : h’*0_{k,n> ; h[ n ; ; g kn and the same expressions with uppercase 

HandG. 

fiCW) . f *{ (CLO)} ; /< c£l > ; p }k : p„{j,k} ; : 

\Long l«f t right arrow 

A*: \la*bda_k ; p k : \au_k ; a\ n : \alpha‘0.{k,n>. 

#? „ : \bwta‘0.-{k,n> ; ©“ : \Thata_j'0; j:.|!w« : M-\I_{H"1>. 
lid/) : \Pi„j(f) ; : \frac{\partial fHApartial x} 

\chi.{[\frac{iHn>,l-\frac{l){aJ]> ; T * : \Xi_j'0. 
fts : \approx ; 
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Table 1: The left and right filter coefficients, //£„ and //* *> for 
the constmction with Divichlet homogeneous boundary conditions 


/(0) = /(l) = 0 . 


Wn 

6.3091928199*01 


1 9.6514114329*02 1.0261005527*01 

2 3.98787 18285e-02 2.2276045032*02 

4 -3.9847422688C-01 -3.8101599048*01 

5 -4.6171223341*01 -3.5419483844*01 

6 -4.3450617898*01 -3.77721 13863e-01 

7 -1.6823324724*01 -4.0626758997e-01 

8 3. 95 12280064c- 02 -2.2139835764*01 

9 7. 43340 18686c- 03 1.5553422358*02 

10 -1.9004116026*02 3.9763630210*02 

k=l 0 -7.1465107925C-01 -6.791 72003S5e-01 

1 1.0487856031*01 7.0564886681*02 

2 -1 0275008552*01 -6.1435155066*02 

4 -2.2668171853*02 1.0664567825*01 

5 -3.4790471529* 01 -6.9557725550*02 

6 -5.4438357904*01 -3. 12669 13055c- 01 

7 2.1013565982*01 -5.5031934226*01 

8 7.1146281575c-02 -3.3164862537*01 

9 8.8701254411*03 2.04 64 047232c- 02 

10 2. 26772204 74c- 02 5.2318055022c- 02 

k=2 0 -2. 2362773 506c-01 -3.4534 868269c-01 

1 5 008265 1220*01 3.9554860545*01 

2 -1.1707940876*01 -5.1663829154*02 

4 7.0092772722*01 -3.4828534215*01 

5 -4 8072499606*02 - 5.6934808461*01 

6 4.1374409814*01 -1.6719040087*01 

7 1.2755670687*01 4.0449494649*01 

8 -6.4$3822028S*02 2.8819962720*01 

9 -4.9723065158*03 -1.4195865812*02 

10 1.2712119110*02 - 3.6292922913*02 


Note; Foe the right edge (t = 1) U»e coeffideoU {Hq w } ut retmted to ibe ^emling fuiKnow 
vrj mud mre luted from right to left The c«ue n=3 cormpoad# to the funciioo 

we have removed mud will therefore aot mppe&r 
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Table 2: The left and right wavelet filter coefficients, G% n and (?£ i0 
for the construction with Dirkhlet homogeneous boundary condi- 
tions /(0) - /(l) = 0. 

<5C — 

k=0 0 2.0 137223334c- 01 2.4372875232c01 

1 5.4229707169e-01 4.5801973839c-01 

2 ■ 7. 0645 62606 7e-01 -3.431 7083287e-01 

4 3.5108313491e-01 7.4312944286c-01 

5 1.5369973568c01 -1.4784113471e-01 

6 -1.0259167455e-0i -1.8898503160e-01 

7 -7.661 542896 le-02 4.538787 1214c- 02 

8 5.5724919877c02 3.3997$23918c02 

9 2.61 85470422c 03 -1. 561321 7459e-03 

10 -6.69*5353810c03 3.991650 1406e- 03 

k»l 0 -6.6848468945*- 14 8.0232416876c- 14 

1 2 8233692882c- 02 3.8 745242 758c- 01 

2 -5.7322775598c 02 -6. 768039234 3c-01 

4 5.9385536874c- 02 -3. 9552887470c 01 

5 -1.1474607116e-01 3.7058026047e-01 

6 -2.9055557195C-0! 2.8173001 140c-01 

7 8.007401031 4c-01 -1.0486l07295c-01 

8 -4.970491 7725C-01 -8.7l04529274c02 

9 -2. 950046684 9e-02 3.44429846 18c- 03 

10 7. 54 204 20519*- 02 8.8056382200c 03 

k=2 0 6.2086522836c- 14 -5. 080389 1 077c 13 

1 -1.480208491 5c- 14 -4.0101204690* 13 

2 6.75O4O52837C-02 1.0243645353c02 

4 -3.6141627882c01 LI 309599230c- 01 

5 7. 864 2373898c- 01 -5.493371Q005c-01 

6 -4.931030677 8e-01 7.6956164180c01 

7 2. 1421713460c- 04 2.8909897939*01 

8 5 6564958243c- 02 -9.033916466?c02 

9 1. 086922 1407e-03 1.23463 124 50c- 02 

10 2.7783077166c-03 3. 1565662736c- 02 

k=3 0 2.5584 258650c- 02 1. 658577 5990e-02 

1 -6.5876?83667c01 -6 8420143939*4)1 

2 -6.8353593454* 01 -6.4585264939c 01 

4 -3.0331 01 2100c-01 2.5993685887c-02 

5 3.7971289157c02 -2.785569441 8c01 

6 6.731 5440450c02 -1.5252127585c 01 

7 1 0944670231c02 S. 7162800751c- 02 

8 7.6074757976* 03 7. 2638444593c 02 

9 -3.8773215636c 04 -2.8584982261*03 

10 9 9126981384c04 7.30 79907303c-03 

Ao<e. For the right edge (x ~ 1) the coefficient* {C$ „} me related to the wavelet e* j,_j 
and me listed from tight to left The esse u— 3 correspond* to the scaling function we h*ve 
removed sod will therefore not appear 
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Table 3: The left and right filter coefficients, and //£*, for 
the construction with Ne^n: aim homogeneous boundary conditions 

/^(O) « / ^(l) = 0, 



n 

. TS. 

w. 

k=0 

0 

7.0710678119*01 

7.0710678119*01 


1 

Q.OOOOOOQOOQe+OO 

0.0000000000*f00 


2 

0.0000000000e+00 

0.0000000000*f00 


4 

3.7878157560*01 

3.3279641360*01 


5 

4.1933439083*01 

3.1763076592*01 


6 

3.9462541967*01 

3.3872841031*01 


7 

1.5279210977*01 

3.6432796801*01 


8 

-3.5885680933*02 

1.9854306805*01 


9 

-6.7511337557*03 

-1.3947818884*02 


10 

1.7259840309*02 

-3.5658770113*02 

k=l 

0 

-6.3999878989*01 

-6.2399400714*01 


1 

2.72561 58164*01 

2.4746950442*01 


2 

2.8255565264*02 

2.1217268985*02 


4 

9.3033115429*02 

3.1423763181*02 


5 

3.9549596240*01 

1.5138263011*01 


6 

5.4799002074*01 

3.4935923414*01 


7 

2,1165351693*01 

5.4440142107*01 


8 

-6.7366756980*02 

3.2169170065*01 


9 

-9.0159055096*03 

-2.0365603789*02 


10 

2.3049919460*02 

-5.2066376096*02 

k=2 

0 

6 3813937293*02 

1.4204323343*01 


1 

6.2433761707*01 

-5.3003493793*01 


2 

7.3965386969*03 

1.7695538529*02 


4 

-6.2421804274*01 

-4.0650580195*01 


5 

-1.3488453082*02 

-4.9029535249*01 


6 

4.3699790930*01 

-9.0099784088*02 


7 

1.4222943082*01 

4.3727314163*01 


8 

-6.9374369289*02 

3.0298274581*01 


9 

-5.5998777298*03 

-1.5509344332*02 


10 

1.4316557613*02 

-3.9650941036* 02 


For the right edge (x = 1) the coefficients {H& K ) are related to the scaling function 
■fj ,. *»*l are lilted from right to left Tbe esse n--3 corrwpood* to the scaling function 

we have removed and will therefore not appear 
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Table 4: The left and right wavelet filter coefficients, 6^ , and G i„, 
for the construction with Neumann homogeneous boundary condi- 
tions /^>(Q ) = /<»>(!) = 0. 



n 

: 

5f7~ 

k--=0 

0 

2. 157828191 4e-01 

2.3494978252*01 


1 

4.9572192648*01 

4.8386498781*01 


2 

7.2238470640*01 

5.2859034242*01 


4 

•39185512760*01 

-6.4278506397*01 


5 

-1.2862537931*01 

4.9429597015*02 


6 

9.6457537590*02 

1.2296776567*01 


7 

6.4796296855*02 

-1.9638518548*02 


8 

•4.7138054627*02 

-1.1782477983*02 


9 

-2.2555618915*03 

7.3127496136*04 


10 

5.7665333650e-03 

1.8695658406*03 

k— 1 

0 

5.2517520677*13 

5.5828490449*14 


1 

2. 4258207630c- 02 

2,0*58047386*01 


2 

6.0376140438*03 

5.9531516311*01 


4 

2.2848095434*01 

5.4090218399*01 


5 

•5.1810004520*01 

4.1031183125*01 


6 

5.6358588265*01 

-3.4180312799*01 


7 

•5.2660919644*01 

1.2529183073*01 


8 

2.8410901284*01 

1.0078064242*01 


9 

2.0215169170*02 

-4.1 780815641*03 


10 

5.1681777360*02 

1.0681616333*02 

CM 

II 

0 

1.9660656197*12 

-2.7202352286*12 


1 

3.0437460103*12 

-9.0382597736*12 


2 

5.5424594564*02 

1.5213229661*02 


4 

2.7483194254*01 

1.2607094728*01 


0 

6.0853476157*01 

-5.5744040614*01 


6 

•1.0850026370*01 

7.6287424658*01 


7 

-6.0447233132*01 

-2.8671529090*01 


8 

4.1282047408*01 

-8.8419425836* 02 


9 

2 1 55405609 1*02 

1.2267369323* 02 


10 

5.5104754193*02 

3 1362559711*02 

II 

0 

1.9941691330*01 

1 8775994937*01 


1 

5.3813012840*01 

6.1793471320*01 


2 

6.8862102186*01 

-6.0431561252*01 


4 

-4.2066297006*01 

•3.7019350637*02 


5 

7.4122879021*02 

-3.8404010349*01 


6 

1.1518683604*01 

-2 0032106281*01 


T 

1 

2.1861286109*02 

1.3094817230*01 


8 

-1.5750923988*02 

1.0651409643*0) 


9 

7.6389913001*04 

-4.3365631441*03 


10 

1 .9529722670*03 

-1.1086787800*02 


Note: For Uw n*hi ed & = eoefficieatt {C| „} are related to the wave!* *' „ . 

and are lifted from rxgbt to left. The tur n^3 coereapoodr to the araling functions we 
have removed and will therefore not appear 
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Tsble 5: Tbe first four moments for the 0 edge scaling function 
of Figure 2 (i.e satisfying Dirichiet homogeneous boundary condi- 
tions). 


X, 

T=*\ 



-1.805l8411015e+00 1.58610287947*41 1 .425773 15923e-01 


1=1 3.4597$M0889e+00 -1. 342791944 17e+ 00 1 66298580543c 13 

1*2 -8. 12722032481*+ 00 - 5.4 1163047176*+ 00 7.89305924956c 01 

1*3 -2.1337583ll07e+01 -1.76739256352e+01 4.70776M5989e+00 


Note: For thus cm* ike monomial x could be expanded as a trw* .ombtMUoc of £ t 
and j. This explains the aero value of the second moment of ,4 } 


Table 6: The first four moments for the 1 edge scaling function 
of Figure 1 (i.e satisfying Dirichiet homogeneous boundary condi- 
tions). 




1=0 j -2.03632572831 *+ 00 2. 11 233673136*41 2 34792863276*41 


1=1 -5. 1 7867977 397c +00 - 1.87776044074c *00 -4 2 1625644274c- 13 
1=2 -1 57657032623e+0l 1.04 701 396775c- 01 1 27839462595c ^00 

1=3 5.24899484894*+Ql -4.53961 860524e+01 1.15515969500*^01 


iVofc; same remarks Table 5 for the value of the tccood moment of , 
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•)The three seating functions for Che loft edge 





0 2 4 6 


b)The thn*e seating functions for tie right edge 



- 6 — 4-2 0 


2. The sir taqe scaling Junctions for the cast \~J. epresented at scale 
j-0 (it. on jO A ex. I for the left edge and j-oo.Oj for the right tdgel Ail 
these scaling functions satisfy Diruhlet homogeneous boundary t'onditwns <j / 
0 
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a)Tbe three scaling functions tof the left edge 





b)The three scaling functions for the nght edge 





FIC». 4 , Sarne a.i in hq 2 jor the Neumann homogentou.* boundary condi- 
tion* 
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a) The lour wavelets lor the left edge 



0 2 4 6 0 2 4 



0246 0246 

b)The four wavelets tor the right edge 



(son . .« 
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